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The need to build a link between the structure of a complex network and the dy- 
namical properties of the corresponding complex system (comprised of multiple low 
dimensional systems) has recently become apparent. Several attempts to tackle this 
problem have been made and all focus on either the controllability or synchronis- 
ability of the network — usually analyzed by way of the master stability function, 
or the graph Laplacian.. We take a different approach. Using the basic tools from 
dynamical systems theory we show that the dynamical stability of a network can eas- 
ily be defined in terms of the eigenvalues of an homologue of the network adjacency 
matrix. This allows us to compute the stability of a network (a quantity derived from 
the eigenspectrum of the adjacency matrix). Numerical experiments show that this 
quantity is very closely related too, and can even be predicted from, the standard 
structural network properties. Following from this we show that the stability of large 
network systems can be understood via an analytic study of the eigenvalues of their 
fixed points — even for a very large number of fixed points. 
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I. INTRODUCTION 



The principal aim of this paper is to ask what does the structure of a complex network 
tell us about the way that it will behave as a dynamical system. We are not the first to 
ask this question. Past attempts have focussed on either conditions to obtain a synchronous 
state, or controlling network dynamics to a target state. Our approach is to take the tools 
of dynamical systems theory and apply these to the network — treating it as a very high 
dimensional dynamical system, with certain useful symmetric properties. We present the 



results of this work in two parts. In Sec. Ill we present a computational study of the 



relationship between various network properties. Sec. |IV| reverts to the basics of dynamical 
systems theory and presents, through several worked examples, direct application of these 
methods to network systems. In the remainder of this introduction and Sec. [IT] we briefly 
review the necessary background concepts. 

Complex networks, and in particular, scale-free networks are widely described as ubiq- 
uitous throughout Nature. While this claim is reasonable, the widely utilized model for 
generating a scale-free network, preferential attachment^, does introduce slight statistical 
bias for finite size realisations^. For a finite network, constructed via preferential attach- 
ment, the connections from the last added nodes are biased towards the hubs, and yet these 
nodes themselves have exceptionally low degree. This leads to finite scale free networks 
constructed with this method exhibiting negative assortativity. To overcome this we have 
proposed a form of altruistic attachment as an alternative^ rather than connect directly to 
the hubs, new nodes are connected to random neighbors of the hubs. By doing this, the 
disassortatvity of preferential attachment is entirely mitigatecP. 

Here, we choose to focus on the effect of this biased assortativity — or disasortativity 
- and ask how does the structural and dynamical behaviour of a strongly assortative (or 
disassortative) network differ from the archetypal preferential attachment model. To quantify 
the structural properties of a network we measure an ensemble of the usual suspects: node 
degree, diameter, assortativity, clustering and robustness (all to be briefly described latter). 
To measure dynamical behaviour we compute the leading eigenvalues of the fixed points of a 
network dynamical system (which is also described more precisely in the following sections). 



Against this battery of statistics we assess four different types of complex networks — two 
standard and two of our own invention: (1) the standard preferential attachment modeP, (2) 
its small-world sibling^, (3) the altruistic attachment model briefly mentioned above, and 
(4) a "skinny" scale-free network designed to be not-small- world 9 . From here on we refer 
to these four types of complex networks as: preferential attachment; small-world; altruistic 
attachment; and, skinny. For each of these four classes of networks we impose a link-exchange 
mechanism (which will also be described below) to incrementally alter the assort at ivity. In 
the next section we introduce the necessary machinery. 



II. NETWORKS AND NUMERICS 



In the following three subsection we introduce the necessary numerical techniques: sta- 



tistical quantification of network structure (Sec. II A), numerical measures of dynamical 



stability (Sec. II B), and methods of link manipulation to modify network assortativity (Sec. 



II C ). In all of what follows we consider unweighted undirected graphs of N nodes represented 
by an A x A symmetric adjacency matrix A such that Ay = 1 indicates the presence of a 
link between node-z and node-j (Ay = otherwise). The main diagonal is zero (An = 0). 



A. Structure 

The measures of network structure which we employ are fairly standard throughout the 
literature, we reiterate their description here briefly and refer the reader tcF' or the relevant 
sources for details. 

• node degree: the average number of links for a node. 

• diameter: the median over all i and j of the shortest path between node-z and node-j. 

• assortativity: the correlation coefficient between the degree of node-i and the degree 
of the neighbors of node-i — computed over all i 

• clustering: the probability that two neighbors of node-i are also neighbours — i.e. the 
prevalence of triangles within the network. 
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• robustness: the tendency for one of the properties of a network (usually diameter) to 
change with targeted or random deletion of nodes — targeted deletion occurs when 
nodes which are deemed to contribute most to that property (high degree nodes,^ for 
the case of diameter) are removed first. 



B. Dynamics 

Let 4>(x) define a dynamical system: x' = 4>(x). In what follows we restrict our attention 
to the one dimensional case and consider a concrete example for the purposes of simulation: 



4>(x) = —x (we will expand on this example in Section IV A). Clearly, this system is globally 
stable with a single stable node^at the origin x = 0. Define the following network dynamical 
system: 

z' = (I + eA)$(z) (1) 

where I is the N x N identity and e > is a small positive coupling strength. The variable 
z consists of N sets of state variable of the system x' = 4>{x) (since, here, x G R, z G R ) 
and $(•) = </>(•) x </>(■) x ... x </>(■) is N copies of the function </>(•). 

If A is full rank, then the solutions of are Zq G {z\(f)(z^) = OVi} where is the 
i — th component of z. In the example we've chosen, this is trivial: z — 0. The stability 
of this single N- dimensional fixed point is determined via the eigenvalues of the matrix 
(I + eA) J$( Zo ) where the Jacobian of $(zo), J<^{z )i ls a diagonal matrix consisting of the 
individual derivatives 4>'(xo) = — 1. For e < 1, (1 + eA) J$( 20 ) ~ — / and these eigenvalues 
are all about —1. Hence, we are interested only in the largest eigenvalue of (/ + eA) J$( Zo ) (or 
equivalently, the smallest eigenvalue of (J + eA)). For a given A we can also determine the 
value of e at which this first eigenvalue changes sign — at this point the system gains a single 
unstable direction, becomes a saddle, and (in the case of <p(x) = —x) becomes unbounded. 



C. Links 



The four networks identified in the introduction exhibit a fairly wide range of assortativity 
- two of them (the skinny network and the small-world network) are highly assortative. 
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FIG. 1. Varying assortativity for preferential attachment networks: The effect of increasing 
and decreasing assortativity via the link exchange mechanism described in the text. Starting with 
the small preferential attachment network (100 nodes) we increase and decrease assortativity with 
50000 applications of the link exchange mechanism to arrive at the assortative and dissassortative 
extrema shown below. At the dissassortative extreme, the network hubs are directly connected to 
the outliers. At the assortative end of the spectrum the nodes are strongly ordered according to their 
degree. The picture is very similar for the altruistic attachment model. The most notable feature 
of small-world and skinny scale free networks is the collapse of the network size with link-exchange. 
We choose very small networks as visualization of larger structures is less informative. 

Nonetheless, how much more assortative (or disaasortative) can a network be? To explore a 
wider range of networks (ordered via assortativity) we implement the following simple link 
exchange mechanism: 

1. pick two nodes i and j at random such that their degrees d(i) and d(j) differ 
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FIG. 2. Varying assortativity for the "skinny" scale free network: The skinny scale free 
network is a rather pathological case — it is a scale free network, but it is not small world. The top 
image depicts the global structure of such a network (of 1000 nodes). Starting from this network 
we increase and decrease assortativity with 50000 applications of the link exchange mechanism to 
arrive at the assortative and dissassortative extrema shown in the lower portions of the figure.. 
At the dissassortative extreme, the network hubs are directly connected to the outliers and the 
network is highly compact. Multiple hubs are also loosely interconnected and are obscured in this 
projection. Conversely, at the assortative end we see community formation based on approximate 
node degree. Note that the network diameter is greatly decreased in either case. 

2. pick neighbors of node-z and node-j, i and j, such that the metrics 

\d(i)-d(j)\ 

and 

\d(j)-d(i)\ 

are minimised 
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FIG. 3. Varying assortativity with a link-exchange mechanism: The six panels depict 
computation of network properties described in the text after repeated application of link-exchange 
to either increase or decrease the assortativity. The four initial networks are: BA preferential 
attachment (red); the WS small- world model (green); "altruistic attachment" (blue); and, the 
"skinny" network (black). Mean and standard deviations from 52 trials are depicted. The horizontal 
axis (for every panel) is the extent to which the link-exchange mechanism was applied to increase 
(postive) or decrease (negative) the assortatitivity. Link-exchange was applied between 100 and 
50000 times (the scale is not linear). 



3. exchange links: disconnect node-i from node-i and instead connect it to node-j, simi- 
larly disconnect node-j from node-j and instead connect it to node-i 

4. repeat 

Note that the individual degrees of each node (and hence the degree distribution) are pre- 
served, but since the degrees of the original nodes differ (d(i) ^ d(j))) and link exchange 
mechanism attempts to produce strong links between nodes with similar degrees, the resul- 
tant network will tend to have a higher assort at ivity. 

Conversely, to decrease assortativity we do the reverse, step [2] is replaced with the follow- 
ing: 

EI. pick neighbors of these nodes i and j such that the metrics 

\d(i)-d(i)\ 

and 

\d(j)-d(f)\ 

are minimised. 

That is, instead of seeking to achieve a good match, we seek to destroy an existing good 
match. 

Over time, these iterative schemes can be applied to significantly increase of decrease 
network assortativity. In the next section we present the results of our computational study 
of these networks. Figure [T] depicts the effect of these exchange mechanisms for the four 
types of networks considered herein. 

III. ASSORTATIVITY AND THE DYNAMICS OF THE NETWORK 

We have performed an extensive set of simulations to fully explore the parameter space 
delineated in the previous section. Figure [3] highlights some of these results. For each of the 
four networks described in Sec. [I] we repeatedly apply the link-exchange mechanisms of Sec. 

8 




FIG. 4. Dynamic stability of the network as a function of network assortativity Each 
network (with the same colour coding as Fig. [3] is parameteerised by its assortativity and maximum 
eigenvalue and plotted here. The right hand panel re-plots results only for the preferential and 
altruistic attachment models. 



II C to increase and decrease the network assortativity. For each network, we then compute 



the range of statistics described in Sec. II A and II B 

From Fig. [3] we can draw several results. The link-exchange mechanism does have the 
desired effect on assortativity and has little effect on scale-free exponent (for networks from 
which this could be reliably estimated). Network degree is invariant — as required by the 
link-exchange mechanism. The link modification seems to have a symmetric effect on both 
path-length and clustering of the small-world and skinny networks. Unlike the preferential 
and altruistic attachment models, both these networks are embeddable, or almost embed- 
dable, in two dimensions (such that nodes that are close according to a Euclidean metric are 
also connected). The link modification scheme disrupts this property and hence decreases 
clustering and path-length — as the networks become more random (exemplified in Fig. [2]). 

Finally, we observe that decreasing assortativity causes a nonlinear increase in the max- 
imum eigenvalue of the resultant network for the scale-free networks — with some notable 
exceptions. More assortative networks are more stable: and the effect is most pronounce 
for the most extremely assortative of disaasortative networks. In Fig. [4] and Fig. [5] we 
probe this relationship further and ask what can one learn about the dynamical structure of 

9 




FIG. 5. Dynamic stability of the network as a function of network mean path-length 
(diameter) Each network (with the same colour coding as Fig. [3] is parameteerised by its mean 
path-length and maximum eigenvalue and plotted here. The right hand panel re-plots results only 
for the preferential and altruistic attachment models. 

the network from examining the structural properties of path-length and assort at ivity. Of 
course, with all else being equal a marginally more assortative network will typically have a 
smaller diameter (path-length) as high degree nodes will tend to be wired together, rather 
than being connected to different regions of the network (the calculation of Fig. [3] illustrates 
this point and the examples in Figs. [T] and [2] demonstrate the opposite effect for larger - 
fairly extreme - changes in assortativity). 

In Fig. [4]we see the effect of assortativity more clearly — increasing assortativity decreases 
the maximum eigenvalue (making the network dynamically more stable). While widely dis- 
tributed, one does observed that for a given value of assort at ivitiy the maximum eigenvalues 
are bounded above by a quantity which decreases linearly with assortativity. At present we 
can offer no good explanation for this linear-in-assortativity upper bound on path length. 
In Fig. [5] we see a similar relationship between mean path length and maximum eigenvalue. 
Taking all three quantities together we find that the maximum eigenvalue can be predicted 
from a combination of assortativity and path-length — in three dimensions the set of points 
depicted here approximates a smooth two-dimensional surface. Figure [6] depicts exactly this 
surface. We build a model (using the modeling procedures we have developed previously^ 
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FIG. 6. Predicting stability Contour plot of a nonlinear model prediction of stability as a 
function of assortativity and path-length overlaid on the data from Fig. [4] and [5] 

to predict stability (as defined above) from just assortativity and path-length. The error in 
these predictions (when applied to new data) has a mean of about 0.025 (~ 12%). 

Finally, in Fig. [7] we compute the effect of changing the network assortativity on the 
robustness of the network. As shown in the upper panel, for all networks, targeted removal 
results in a linear increase in network size (diameter) as the removal fraction increases — up 
to the point where the network fractures. Random removal has very little effect on network 
size. This effect is robust across all levels of initial assortativity. Hence, the robust and fragile 
property of these networks is preserved. Conversely, targeted node removal has a very marked 
effect on network assortativity — by removing a small fraction of nodes assortativity is very 
quickly attenuated. Targetted link removal attenuates assortativity (making disassortative 
networks more disassortative and assortative networks more assortative) for all three classes 
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FIG. 7. Network robustness Results for the preferential attachment model are shown - the 
altruistic attachment model produces qualitatively identical results, the other networks are more 
idiosyncratic. Again, we apply the link-rewiring scheme to change the assortativity of the network 
over a range of values — the same range of values as depicted in Fig. [3| from highly assortative 
to disassortative (as the axis labelled "assort." increases). For each such network we then perform 
targeted (the taper plot) and random (the translucent surface) node removal. Computed values of 
network diameter (mean minimum path-length) are shown on the upper panel and assortativity on 
the lower. 
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of scale-free networks. Again, this is a curious observation and worthy of closer examination 
at a future time. The effect is least pronounced in the skinny network and evident in the 
non-scale-free small world model only at the smallest fraction of link removals. With targeted 
removal of more links, small-world networks become less assortativity — as the structure 
inherent in the network is disrupted. Across all types of networks, targeted link removal 
leads to an increase in clustering (since clustering indicates redundant links, and targeted 
link removal avoids such links). 

The details of this computational study highlight several interesting connections among 
the various topological properties of complex networks — and connections to network dy- 
namics. The basic message from this work is that structure (and perhaps especially assor- 



tativity) matters. In Sec. IV we turn to looking more closely at exactly what this means for 
the dynamics of these systems. 



IV. ANALYSIS OF STABILITY OF NETWORKS 



To define the measure of stability we used throughout the previous section we introduced a 



specific network dynamical system in Sec. II B We now revisit and expand on that definition 



here. In the subsequent subsections we work through several instructive examples. 

As before, x' = <fi{x) denotes the dynamics on each node and A be the network adjacency 
matrix (hence Ay = 1 iff node-z and node-j are linked and = otherwise). Previously, 
we assumed that x G R, we now relax that assumption and allow the dynamics on the node 
to be a system of k differential equations- 7 x E R fc . Construct $(■) = </>(■) x </>(■) x . . . x </>(•) 
from N copies of <j). 

For the sake of concreteness we assume that the nodal coupling is positive and bidirec- 
tional. Moreover, we assume that coupling occurs only amoung the first components of x - 
the only reason for this assumption is to simplify the construction of the nk-bj-nk coupling 
matrix B, below. Although the matrix B is induced from A under the restrictions on cou- 
pling just described, the analysis that follows will work equally well for an arbitrary matrix 
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B. Let 




A 



o 



if both i and j are integer multiples of k 



otherwise 



(2) 



For one dimensional nodal dynamics, B = A, otherwise it is a sparse (sparser than A) matrix 
indicating coupling between the first component of the individual sub-systems, (k — 1) of k 
rows and columns are all zero. 

The global complex system dynamics can now be described by 



where (again) e > is the coupling strength and the dynamics on the n-th node are described 



The natural first step towards understanding a dynamical system such as (J3J) is to find 
the fixed points, linearise them and then use this to characterize their stability. If zq is a 
fixed point we have (/ + eB)^(zo) = 0, i.e. $(z ) G {u\(I + eB)u = 0} the nullspace of 
(I + eB). If (/ + eB) is rank deficient, this then leads to a vector sub-space of (/ + eB) for 
which any solution of §{z) on that subspace is a fixed point. Such a non-trivial nullspace 
could (somehow) be factored out, so let us assume that the only solution is the trivial one 
and that (/ + eB)Q(zo) = implies $(^o) = 0. Indeed, this is not a particularly onerous 
restriction, since a non-trivial nullspace implies nodes in A (and B) that can be expressed 
as a linear combination of one another and could be trivially condensed to a single node. 

Hence, without loss of generality, the fixed points of pi) are the fixed points of x' = 4>{x). 



That is, zq is a fixed point iff 4>(zq ,Zq , . . . , Zq ) = for n = 1, . . . N. We will 



come back to this later, but notice that this introduces quite a bit of symmetry — if 4>(x) = 
has m solutions then their are m N fixed points in ([3j situated at points in an iV-dimensional 
hyper-lattice. Moreover, the location of the fixed points is independent of B (and A) - 
provided that no nodes in the network are redundant (i.e. A has only the trivial nullspace). 

Suppose that <p(x) = has m solutions: a±, a%, . . . a m and call the set of these solutions 
A = {ai, a 2 . . . a m }. Then z is a fixed point of (|3]) iff zq e A x A x . . . A. Hence, which fixed 
point we are interested can be identified with an iV-tuple of integers from 1 to m. Let £(t) 



z = 



(I + eB)<S>(z) 



(3) 
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denote the value of the i-th component of that iV-tuple — hence the fixed point we are inter- 
ested in from the t-th subsystem is am) . Moreover, we can write z = [a^ a^ 2 ) ■ ■ ■ a^jv)] ■ 
Let J^m) = bi and then the stability of the fixed point z is characterized by the eigenvalues 
and eigenvectors of the matrix 



6/(i) 














6/(2) 




























6/(jv) 



(4) 



Computationally this is all perfectly do-able, and the symmetry we alluded to earlier actually 
makes it quite simple for many systems of interest. 



A. Single wells: 1 fixed point 



In Sec. II B we considered the simplest possible non-trivial dynamical system. We do it 
again here. Let <fi(x) = —x. Hence, for e = the dynamics on each node consists of a one 
dimensional stable node at xq = 0. This one dimensional system means that A = B. As we 
increase e « 0, the matrix (I+eA) > has only positive eigenvalues. Hence the eigenvalues of 
Q are all negative — the system has a single stable node at the origin zq = 0. Eventually, 
for e large enough the matrix ^ will gain first one and then increasingly more positive 
eigenvalues (two positive eigenvalues could emerge simultaneously given certain degenerate 
- and symmetric — arrangements of node in the network). But, as far as bounded dynamics 
are concerned, one unstable direction is sufficient: at this point the node becomes a saddle 
and all trajectories diverge. 

At first this may seem a little counter intuitive, and hence is interesting. Each node is 
very stable and yet by coupling these nodes together the system dynamics finds a way to 
climb out of this one dimensional potential well — diverging to ±oo. Moreover, despite the 
universally positive coupling, some nodes diverge to oo while other go to — oo., Figure [8] 
depicts typical trajectories, confirming this behaviour and validating our choice to employ 



the largest eigenvalue of Q as our measure of stability in Sec. Ill 
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FIG. 8. Network robustness Sample trajectories (the trajectories for the 10 highest connected 
nodes are shown, all others are less divergent) for A a preferential attachment scale-free network, 
0(x) = -x and N = 1000 nodes. Top-left is with e = 0.092 and top-right is e = 0.093. The 
eigenspectrum of I + eA is shown in the lower plot (the eigenvalues of ^ have the opposite sign 
because <f>'(x) < 0). At e = 0.092 all eigenvalues are negative, at e = 0.093 one becomes positive. 



Finding the critical value e c such that the system's fixed point gains one unstable direction 
is also straightforward. Rewrite (|4j) in this case as ((— eB) — I) since bg^ = (j)'(xo) = — 1. 
The eigenvalues of this matrix are the solution of det ((—eB) — (1 + A) I) = 0. Denote by 
Xb the eigenvalues of B — which satisfy det (B — XbI)- Hence A = —1 — e\s and the 
eigenvalues of our single potential well network system will all be negative if 1 + eA# > 0. 
Hence the critical value e c for the node-saddle bifurcation is e c A# > — 1 for all the eigenvalues 
of B. Since e > this is always true of the positive eigenvalues. Hence we require e < —j^. 
This is first violated by the negative eigenvalue of largest magnitude, i.e. \B,mm- Hence, 
the first bifurcation from node to saddle occurs at e c = — , where As m i n is the smallest 
eigenvalue of B. Subsequent transitions could be similarly computed, at least in principle. 
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FIG. 9. Network robustness Sample trajectories for A a preferential attachment scale-free 
network (the same network as that used in Fig. [8]), 4>{x) = x(9 — x 2 ) and N = 1000 nodes. Top-left 
is with e = 0.092; top-right is e = 0.093; bottom left is e = 0.094; and, bottom right is e = 0.097. 
At e = 0.99 the system becomes unbounded. The eigenspectrum is the same as that in Fig. [8j 

B. Double wells: 3 n fixed points 

The next situation is a little more interesting. For the time being we will stick to one 
dimensional dynamics, but we now introduce slightly more complex dynamics on the nodes. 
Let 4>(x) = x(9 — x 2 ). The system has three fixed points xq = ±3,0. In the nomenclature 
introduced above CLi — — 3, &2 = 0, and = 3. The values of <\>' are sufficient to characterize 
the stability of the fixed points b\ = 0'( a i) = — 18, b 2 = 0'( a 2) — 9, and 63 = 0'( a 3) = — 18 ; 
I.e. two stable and one unstable nodes. 

Again, for a small e > the positive definiteness of (/ + eA) ensures that the fixed points 
behave independently of one another. Hence the fixed points of ^ occur at zq where each 
component of zq is —3, or 3. These are the vertices of a hypercube in H N together with 
the mid-points of each edge and face and the centre (at the origin). There is one unstable 
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node (corresponding to the fixed point at 0) and 2 N stable nodes (each vertex of the cube). 
All other fixed points are saddles. As before, once e increases the stable fixed points at the 
vertices of the cube gain unstable directions and become saddles — computational simulation 
verifies that soon after this happens the system is no longer bounded. Because 4>'(xq) is a 
fixed constant at the stable nodes, this critical value is the same as in Sec. 



IV A 



Figure [9] confirms this analysis. Moreover, further computational probing shows, in Fig. [9j 
that the system develops first a stable periodic orbit around the origin and then a stable focus 
at — that is, the fixed point which is initially unstable now becomes stable. Note that the 
fixed point depicted in Fig. [9] is not the original in R^, but a point z = [a^i) a^ 2 ) • • • a^(jv)] 
such that there exists a unique i such that e^) = and for all other j 7^ i, e^j) = ±3 - 
that is, this is a point on one edge of the hypercube of fixed points described above. 

Finally, the more complex dynamics are a result of the correct orientation of stable and 
unstable manifold of saddle points and hence depend sensitively on the choice of the matrix 
A — not just on the type of the network. 

C. Oscillators 

So, what happens if you start which oscillatory node dynamics? Can probably subsume 
this into the next example. 

D. Networked hyper-chaos: higher dimensional nodal dynamics 



A lot of the notation we introduced at the beginning of Sec. IV is unnecessary if (f)(x) is 



a scalar function. We now attempt to make better use of that notation. Let x = (xi,X2,x 3 ) 
and 

/ . . \ 



x 2 
V s J 



(5) 



—x 2 - x 3 
X\ + ax 2 
yb + x 3 (x 1 -c) J 

where b = 2, c = 4 and a is a bifurcation parameter. This is, of course, the Rossler system 
where for a = 0.398 the system exhibits broad band chaos and for smaller values of a 
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FIG. 10. Network of chaotic oscillators: Sample trajectories for A a preferential attachment 
scale- free network (the same network as that used in Fig. [8]), <p(x) given by ^ and N = 1000 
nodes. Various values of the bifurcation parameter a and the coupling strength e are illustrated. 
Row-by-row, the uncoupled dynamics of the Rossler system at that value of a exhibits: (1) limit 
cycle; (2) period 2; (3) period 4; (4) "four-band" chaos; (5) period 6; (6) "broad-band" chaos; (7) 
period 5; and (8) period 3. An expanded view dfithe system for a = 0.398 and e = 0.02 is shown 



in Fig. 11 




FIG. 11. Driven chaotic Rossler: Sample trajectory for A a preferential attachment scale-free 
network (the same network as that used in Fig. [8]), <p(x) given by j}, N = 1000 nodes, a = 0.398 
and e = 0.02. We plot the dynamics on the largest hub for 10000 time steps after a transient of 
10000 time step (sample time 0.1). 

exhibits periodic dynamics and a period-doubling bifurcation via the archetypal stretching 
and folding mechanism. 

Although an analytical treatment of N = 1000 coupled Rossler systems may be a little 
difficult, we can, at least, study the effect of such a coupling computationally and relate it to 
the dynamics we have already observed for simpler systems. Figures 10 and 11 summarise 
those calculations. Note that here e < e c — the interesting dynamics highlighted here occur 
prior to any change in the stability of the fixed points. 
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V. CONCLUSIONS 



It is not sufficient to collectivise all scale free networks and treat them equally. While 
some properties — such as the so-called "robust-yet-fragile" character of these networks is 
preserved upon manipulation of assortativity. Other properties are sensitively dependent 
on the local connection between nodes. In particular, more assortative networks are more 
dynamically stable. By merely changing the assortativity on a network one can progress from 
a regime of stable equilibrium dynamics to global instability. Nonetheless, we do find that, 
to a very large degree, stability of a network can be predicted from just network diameter 
and path-length — two global properties of connectivity that are, in turn, strongly affected 
by assortativity. 

Treating the networks as dynamical systems we are able to make use some basic results 
from dynamical systems theory to obtain some level of understanding of the evolution of 
complex dynamics in these systems. We see that even for a fairly large dimensional systems, 
if the dynamics on then individual nodes are identical then the system can be understood in 
terms of those fixed points. Moreover, the initial transition from a stable node to a dynamical 
saddle does not depend on (f)'(x ), but can be determined directly from only the smallest 
eigenvalue of A. 

In computational simulations of large network systems we observe interesting effects of 
the compounding dynamics. Single well-potential functions can be forced to "climb" away 
from the stable equilibrium. Unstable equilibira in double wells are stabilised, and for chaotic 
systems the dynamical coupling appears to intensify the complexity (and possible the chaotic- 
ness) of the nodal dynamics. 

Finally, we computationally demonstrate a strong dependence of the system dynamics 
(characterized by this smallest eigenvalue of A) on the median path-length and assortative 
of the network. Moreover, we have uncovered an apparent linear-in-assrotativity upper 
bound on the dynamical stability of these network systems. 
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